Geometric phases and Bloch sphere constructions for SU(N), with a complete 

description of SU(4) 
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A two-sphere ("Bloch" or "Poincare") is famiUar for describing the dynamics of a spin- 1/2 particle 
or light polarization. Analogous objects are derived for unitary groups larger than SU(2) through an 
iterative procedure that constructs evolution operators for higher-dimensional SU in terms of lower- 
dimensional ones. We focus, in particular, on the SU(4) of two qubits which describes all possible 
logic gates in quantum computation. For a general Hamiltonian of SU(4) with 15 parameters, 
and for Hamiltonians of its various sub-groups so that fewer parameters suffice, we derive Bloch-like 
rotation of unit vectors analogous to the one familiar for a single spin in a magnetic field. The unitary 
evolution of a quantal spin pair is thereby expressed as rotations of real vectors. Correspondingly, 
the manifolds involved are Bloch two-spheres along with higher dimensional manifolds such as a four- 
sphere for the SO (5) sub-group and an eight-dimensional Grassmannian manifold for the general 
SU(4). This latter may also be viewed as two, mutually orthogonal, real six-dimensional unit vectors 
moving on a five-sphere with an additional phase constraint. 

PACS numbers: 03.67.-a, 02.20.Qs, 03.65.Vf, 03.65.Fd, 02.40.Yy 
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INTRODUCTION: THE BLOCH SPHERE 
AND ITS EXTENSION 



In the study of the dynamics of a spin-1/2 particle, a 
visual metaphor that has played a powerful role is that 
of the "Bloch sphere" Pure states of the system are 
represented by the tip of a vector from the origin to the 
surface of such a unit sphere S^. In the field of nuclear 
magnetic resonance (nmr) Q and elsewhere, transfor- 
mations between states are then viewed as rotations of 
that vector, described by the Bloch equation of motion, 
m — —2B X m, for a magnetic moment in a magnetic 
field B. Thus, various sequences of nmr manipulations 
can be pictured in a nice geometrical way as successive ro- 
tations, and this has now become central to our intuition 
of spin dynamics. The relevant group of unitary trans- 
formations is SU(2), a rank-one, three-parameter group 
that is the double covering group of the three-dimensional 
rotation group SO (3) Q. The three operators of angu- 
lar momentum, ( J^,, Jy, J^), are the generators of these 
groups. A canonical set of parameters of S0(3) are the 
Euler angles. Integer values j = 0, 1, . . . provide various 
(2j + l)-dimensional representations, while for SU(2), the 
half-odd integers occur as well. The two coordinates on 
S^, together with a phase, provide the three parameters 
describing the full state. 

This latter phase is often not accessible as, for instance, 
when dealing with the density matrix p. Mixed states 
also are naturally accommodated in this picture. They 
are represented by points inside the sphere so that the 
vector is of length less than unity. Correspondingly, Tr 

< Tr p, which constitutes a definition of a mixed state 
States of light polarization, also a two- valued object, 
map onto the same mathematics and geometry through 
the "Poincare" sphere 0]. 

It would be of interest to have analogous geometri- 
cal pictures for multiple spins, especially in today's fields 



of quantum computation, cryptography, and teleporta- 
tion, because the fundamental elements of these subjects 
are built up of a few qubits @. Thus, all logic gates 
for quantum computation can be built up from qubit 
pairs, while telcporting one qubit state requires an en- 
tangled pair held by the sender and receiver, for a total 
of three qubits. With SU(2p) being the relevant group for 
p qubits, this calls for a similar geometrical description of 
higher SU(7V). In this paper, we develop such a picture, 
through an easily accessible procedure which iteratively 
descends from N to N — n, with n < N, in a manner 
that closely follows the description of SU(2). 

Our procedure also applies when N is odd, a situation 
that does not arise with qubits but elsewhere widely in 
physics (for example, qutrits ?|, neutrino oscillations [8], 
the quark model and quantum chromodynamics (QCD), 
etc.) The N x N matrices of Hamiltonians and evolution 
operators are viewed as built up of 2 x 2 block matrices 
through this N = {N — n) + n decomposition, the block 
matrices then described in terms of the Pauli spinors of 
SU(2). Each step of this iterative reduction introduces 
an analog of the Bloch sphere, albeit of higher dimension 
and more complex structure, and constructs the effective 
Hamiltonians of dimension (N — n) and n for the next 
step. In this manner, using no more than the operations 
familiar from the SU(2) case, the full construction for 
SU(7V) is achieved. 

The philosophy behind such a construction may be 
seen as generalizing Schwinger's philosophy for represen- 
tations of SU(2) or S0(3), where higher j-representations 
are constructed from those of the fundamental, j — 1/2 
[11. We now do the analogous step of using SU(2) as the 
template for solving larger SU(7V). In particular, for the 
important case of SU(4) for two qubits, we give a com- 
plete description of the manifolds and phases involved 
and analytical expressions for them. Note again, as with 
light polarization and spin-1/2, that the mathematics of 
A^-level systems in quantum optics, atomic and molecular 
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physics, and elsewhere, is the same as that we describe 
in the language of multiple qubits. This provides an even 
wider context for our results. 

The arrangement of this paper is as follows. Section II 
describes the basic iterative decomposition of the evolu- 
tion operator for SU(A'^), mimicking the familiar proce- 
dure for spin-1/2. With = 4, and n = 2, Section III 
specializes the results to SU(4), the case of two qubits, 
when all the manipulations involved are in terms of Pauli 
spinors. It also applies these results to Hamiltonians in- 
volving a restricted set of operators of the full group. An 
interesting one is SO (5), which can be described by a 
5x5 antisymmetric matrix that is the analog of the 3x3 
antisymmetric one for the magnetic field B in the Bloch 
equation. Section IV then considers Hamiltonians re- 
quiring the full SU(4) group for their description. Linear 
equations, analogous to the Bloch equation, are derived 
in terms of vectors m, five- and six-dimensional vectors, 
respectively, for the S0(5) and full SU(4) cases. The 
latter also correspond to so-called "Pliicker coordinates" 
[lOj which are also presented. Appendix A deals with the 
generalization to non-Hermitian Hamiltonians, and Ap- 
pendix B presents the isomorphism between SU(4) and 
the groups Spin(6) and S0(6) which we exploit. 



II. ITERATIVE CONSTRUCTION OF 
EVOLUTION OPERATOR IN TV DIMENSIONS 

We wish to obtain the evolution operator JJ^^\t) for 
the iV-dimensional time-dependent Hamiltonian H*^^^: 



H(^-")(t) V{t) 
VHt) H(")(<) 



(1) 



We have blocked the Hamiltonian into {N — n)- and n- 
dimensional blocks, the diagonal blocks being square ma- 
trices while the off-diagonal V is {N — n) x n and is 
n X {N — n). Although our discussion is for Hermitian 
H'^^^ the procedure can also apply more generally, in 
which case the off-diagonal blocks will not be simply re- 
lated as adjoints (see Appendix A). We will also assume 
H*^^^ to be traceless, again a restriction that can be eas- 
ily relaxed, the time integral of the trace becoming an 
overall phase of U*^^^. 

To solve the evolution equation, with an over-dot de- 
noting derivative with respect to time. 



^U(^) {t) = H(^) (i)U(^) (t), U(^) (0) = I, (2) 

we similarly block the unitary matrix, writing it also as 
a product of three factors, the first two further grouped 
as Ui and the second, U2, block-diagonal in form: 



uW(t) 

Ui 



U 



2 — 



U1U2, 

jiN- 

ot 

U(^-»)(i) 
ot U(")(t) 



z(i) 
l(«) 



(3) 



where A± are matrix generalizations of the Pauli spin 
step- up/down a±, and z and are rectangular matrices 
of complex parameters. 

The above structure, with Ui having blocks of zero 
in the lower and upper off-diagonal blocks of its matrix 
factors, is crucial in our method. For the case of spin- 
1/2 and SU(2), the form of a product of three factors, 
each an exponentiation of one of the Pauli spinors, is 
well known Their Cartesian form, with Euler angles 
in the exponents, is the familiar choice but we choose in- 
stead the triplet, (cr±, Cz), when the first two factors have 
zero off-diagonal entries. This introduces complex z and 
in place of the Euler angles, and makes the individ- 
ual factors in Eq. ^ not separately unitary although our 
construction ensures unitarity of the full U'^^'(<). Fur- 
ther, for non-Hermitian H when U is non-unitary, our 
construction still applies. The specific structure of an 
upper and lower triangular matrix and a diagonal one 
proves fruitful, giving simpler equations for z and w''', 
which will have at most quadratic nonlinearity in these 
parameters and not more complicated trigonometric de- 
pendences as with the Euler angle decomposition pTl[T3|. 
They also yield more naturally to a geometrical picture 
of the manifolds they describe. 

A remark about notation. We will use the symbol tilde 
when the corresponding Hamiltonians or evolution oper- 
ators may not be Hermitian or unitary, respectively. Uni- 
tarity of the full U'^''(t) leads to relations between z and 
which would otherwise be independent for evolution 
under a non-Hermitian Hamiltonian (see Appendix A), 



z = -W72 = -7i w, 
7i = U(^-")u(^-")t ^ l(^-»)+zzt, 

72^1 = U("^U(")^ = +z^z)-i. (4) 

With U = U1U2, Eq. ^ formally reduces to the evo- 
lution of U2 alone with an effective Hamiltonian [l^, , 

i\J2 = HoffUa, Heff = Ui 'HUi - iUi ^Ui. (5) 

A key element of our construction lies in this effective 
Hamiltonian and corresponding evolution for the reduced 
problem. Since U2 and this equation are block diag- 
onal, the off-diagonal blocks in Hcs on the right-hand 
side must vanish. This condition leads to the defining 
equation for z. 



H(^-")z + V - z(Vtz + H(")). 



(6) 
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For SU(2), when A'' = 2, n = 1, all the matrices above 
reduce to single numbers and Eq. ([S]) is a Riccati equa- 
tion for the complex z. More generally, it is a matrix Ric- 
cati equation [14||, and its solutions are involved in the 
subsequent construction. With the off-diagonal blocks 
of Eq. accounted for, the diagonal ones defining the 
Hamiltonians for the {N — n) and n problems remain, 
and are given by (H^^-") - zVt) and (H(") -t- Vtz), re- 
spectively. Although the overall trace is preserved in our 
construction and remains zero, these individual Hamilto- 
nians are neither traceless nor Hermitian. The equations 
for z need to be solved numerically in general but form a 
smaller set than the N'^ elements in the original Eq. ([2]). 

To set up the process for iteration, the above individual 
Hamiltonians in (N — n)- and n-dimensional subspaces 
must be rendered Hermitian and traceless. The latter is 
easily achieved, by subtracting Tr (H^^^"^ — zV^) and 
Tr (H^") + V^z) from them. These traces being equal 
and opposite, this translates into the introduction of a 
phase, the integral of the trace, in U'^^-' , representing a 
relative phase between the two subspaces. 

There are alternative methods for rendering the Hamil- 
tonians Hermitian, the most accessible one being through 

The first part of this equation is the observation that 
UjUi is block diagonal. This suggests the second part of 
the equation, namely, the definition of an inverse through 
two "Hermitian square-root" matrices g^. Together, they 
serve as a gauge factor to unitarize according to 

With that, the second factor, U2, in Eq. ((3]) is also uni- 
tarized. 

After some algebra, the explicitly Hermitian forms of 
the two diagonal block Hamiltonians of dimension (N—n) 
and n are 

H(^-") = ^ [Isi"'' +1 (gr' (H^^-^- -Vt) gi+hc) , 

H(") = ^ [^g^"'' §2]+^ (g2"'(H(")+ ztv) g2+hc) ,(10) 

with commutator brackets in the first term, and he in the 
second term denoting the Hermitian conjugate of the pre- 
ceding expression. Again, the trace of each Hamiltonian 
in Eq. (fTU|) can be subtracted to render them traceless; 
as clear by inspection, this is the same trace discussed 



just above. These Hamiltonians in Eq. (fTO|) can now be 
treated further as SU(A^ — n) and SU(n) problems. 

The 7 matrices in Eq. (j4|) are Hermitian with non- 
negative eigenvalues because of their origin from U|Ui. 
This permits their decomposition into g as shown in 
Eq. ([7]). The g matrices and their inverses in Eq. ([7|)- 
Eq. (jlOp . are square roots of them, and because any 
power, including fractional ones, are Hermitian term by 
term in a formal power-series expansion, we can choose 
g also as Hermitian. The use of identities such as 

zt7,P = 7/zt, 7iPz = z7/, (11) 

serves to express all g in terms of the linearly indepen- 
dent set of matrices of dimension {N ~n) or n, whichever 
is smaller. With n = 2, this means that all the algebra 
of calculating such square-root matrices and the subse- 
quent evaluation of the effective Hamiltonian in Eq. (fTO|) 
reduces to manipulation of Pauli matrices. 

A count of the parameters is instructive. The original 
SU(7V) evolution involves (iV^ — 1) elements and, there- 
fore, grows quadratically with N. These are divided in 
the above construction into the 2n{N — n) elements in 
z, which for small n grows only linearly with N. The 
rest are contained in the elements of the SU(Af — n) and 
SU(n) and the single phase between those two subspaces. 
Our construction of higher SU(7V) evolution in terms of 
smaller ones, with the template in Eq. ^ of three factors 
as in SU(2), resembles the Schwinger scheme of generat- 
ing higher j representations of SU(2) or SO (3) from the 
fundamental one of j = 1/2 Q. Whereas that scheme 
was for higher representations but of the same group, 
SU(2), our procedure extends in the direction of larger 
groups SU(7V). 

In mathematical language of base manifolds and fiber 
bundles [15|], the SU(2) and its Bloch sphere are seen 
as the bundle [SU(2)/U(1)] x U(l), the former the two- 
sphere base and the latter U(l) phase the fiber. Like- 
wise, our construction is in terms of the base manifold 
[SU(7V)/(SU(iV - n) X SU(n) x U(l))] and the fiber 
(SU(7V - n) X SU(n) x U(l)). For SU(2), there is a 
single complex z that defines the base manifold. The 
Bloch sphere of a unit three-dimensional vector m cor- 
responding to z is then constructed by inverse stereo- 
graphic projection from to S^. Similar structures of a 
TO associated with the larger z will be considered in the 
next sections. 



III. THE CASE OF SU(4), WITH 
APPLICATION TO ITS SUB-GROUPS 

An important case is of A^ = 4. Four-level systems are 
commonly considered in quantum optics and molecular 
systems and, of course, in today's quantum computation 
where they describe two qubits [1]. Since all logic gates 
can be built up from such qubit pairs, the study of the 
evolution operator for such = 4 problems is of current 
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interest. As a combined description of spin and isospin, 
SU(4) also has central importance in the study of nuclei 
and particles The group also occurs in the descrip- 
tion of unusual magnetic phases of / electron states in 
CeBg |17J . Both choices n — 1 , 2 in the general procedure 
of Section II lead to interesting decompositions, with the 
latter the more natural for qubit applications. We now 
turn to this case. 

In physics terms, a 4-level Hamiltonian has three real 
parameters along the diagonal to fix the energy positions 
of the levels. (One overall element, represented by the 
trace, can be subsumed as an uninteresting definition of 
the zero energy reference level, leading also to an irrel- 
evant overall phase in the evolution operator.) In addi- 
tion, six off-diagonal couplings, which are complex, make 
for a total of 15 parameters to describe the full Hamilto- 
nian. Symmetries often reduce this number so that the 
Hamiltonian involves only a smaller number as a closed 
sub-algebra. For two identical qubits, there are indeed 
such symmetries which reduce the number of indepen- 
dent energies and couplings of a 4-level system. 

With N — A,n = 2, all the matrices involved in the 
previous section can be rendered in terms of Pauli spinors 
and the unit 2x2 matrix. A general Hamiltonian of 
SU(4) has 15 independent operators and time-dependent 
parameters multiplying them. A standard, explicit ren- 
dering of the 15 4 X 4 matrices is given in jisl . [T9| . z 
comprises four complex quantities, (24, z^), and the ma- 
trix Riccati equation reduces to coupled first-order equa- 
tions in them with quadratic nonlinearity. Deferring this 
general case to the next section, we consider first the 
smaller sets of operators of various sub-groups of SU(4). 

su(2) X su(2) sub-algebra: Consider first a Hamilto- 
nian consisting of only six of the 15 operators. Since our 
construction is representation independent, in a suitable 
representation, the six may be viewed as two indepen- 
dent, mutually commuting, triplets that obey su(2) al- 
gebra. Clearly, each then may be expected to have its 
own geometrical description in terms of a Bloch sphere 
and phase. In our above, general formulation, this re- 
sult is realized as follows. Thus, consider two indepen- 
dent magnetic moments, characterized by the standard 
Pauli matrices tr, in time-varying magnetic fields A{t) 
and B{t) which may also be independent, with Hamil- 
tonian H = (T^^' • A + (t'^) • B. Using a standard set of 
4x4 matrices ^8j] to cast this Hamiltonian in the form of 
Eq. (H]), we have V = {A^-iAy)l and Hf^^^) ^ g.B±A,l. 
The z in Eq. ^ also reduces, as with V, to a unit opera- 
tor with a single complex coefficient Z4 obeying a Riccati 
equation in Eq. ([5]). The gamma matrices in Eq. ^ are 
also proportional to the unit operator, thus simplifying 
Eq. pO)) . the g dropping out. As a result, the Hermitian 
matrices in the block-diagonal effective Hamiltonian take 
the form of the same a-B plus/minus a term proportional 
to a unit matrix. The first term is viewed as for a single 
spin with a Bloch sphere and a phase, the second repre- 
sents a phase between the two 2x2 spaces. The complex 
Z4 can again be inverse stereographically projected into 



another two-sphere as in the Bloch construction. We ar- 
rive, therefore, at the same initial expectation, that a 
simultaneous viewing in terms of two Bloch vectors in 
individual two-spheres, along with their fibers, provides 
the geometrical picture for all such qubit-pair systems. 
A specific physical example occurs in the construction of 
optimal quantum NOT operations [20^]. 

su(2) X su(2) X u(l) sub-algebra: Another sub- 
algebra, involving seven of the 15 operators, has been 
considered before It has the symmetry of SU(2) 

X SU(2) X U(l). In a suitable representation, such a 
Hamiltonian can be cast as a diagonal form in Eq. ([T|) 
plus a term which is proportional to the unit operator in 
both diagonal blocks but with equal and opposite sign. 
Such an operator commutes with all the other six, them- 
selves comprised of two mutually commuting triplets of 
4x4 matrices 118]. With V = 0, z also vanishes and 
we reduce trivially to the two independent SU(2) and a 
phase between the two spaces, together accounting for 
the 7 parameters of this problem. An example is pro- 
vided by the CNOT gate constructed with two Josephson 
junctions [l^]- Many such sets of seven operators, one of 
which commutes with all the remaining six, have been 
identified through a general procedure in footnote 11 of 

so(5) sub-algebra: Proceeding further to other sub- 
groups, a non-trivial example is provided by a -ff that 
involves ten operators satisfying an so(5) sub-algebra of 
su(4). Again, there are many such sets of ten opera- 
tors/matrices which close under commutation within the 
full set of 15 as noted in footnote 11 of 21]. As a physical 
example, a four-level system of two symmetric pairs, as 
naturally so with two identical qubits, has only two real 
parameters along the diagonal in its H. Selection rules 
often restrict the off-diagonal coupling between the levels 
from six to four, thus introducing four complex, or eight 
real, parameters. The net result of such symmetric four- 
level systems is a ten-parameter problem [2^. Such H 
fall into this so(5) sub-algebra. The corresponding group 
is the so-called spin group Spin(5) which is the double- 
covering group of SO (5), the group of five-dimensional 
rotations, much as Spin(3), isomorphic to SU(2), is the 
covering group of S0(3) [2J]. All such Spin(5) or S0(5) 
will themselves have a Spin(4) or SO (4) sub-group, which 
in turn has the two mutually commuting SU(2) or SO (3) 
discussed above so that the ten matrices can be conve- 
niently viewed as two sets of commuting triplets plus four 
more which transform like a four-dimensional vector un- 
der S0(4). For completeness here in this paper, we briefiy 
summarize results on this so(5) sub-algebra that were 
published elsewhere [13]; see also psj. 

In a convenient representation that uses Pauli matrices 
for two spins [11], we have H{t) — F2iai^^ - F^iuf''' -\- 
F32<Jx^ -F^^a^P +F^ia^^^ af^ -F^ia^y^ , where the ten 
arbitrarily time-dependent coefficients F^i, (t) form a 5 x 5 
antisymmetric real matrix. (We will use fJ.ji' — 1 — 5 and 
i,j,k= 1 — 3 and summation over repeated indices.) Sev- 
eral quantum optics and multiphoton problems of four 
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levels driven by time-dependent electric fields have such 
a Hamiltonian. It has also been considered extensively 
in coherent population transfer in many molecular and 
solid state systems [111- Casting this Hamiltonian in the 
form of Eq. ([T]) , we have 



a U(l) phase, parallel the discussion of the general SU(4) 
in the next section, and will be presented elsewhere [26j . 
Therefore, we omit details except to note that setting 
Z4 — —izi and 2:3 = —iz2 in Section IV reduces to such a 
SU(3) symmetry. 



jj(i,2) ^ ^^p^^ _ ie,^.^F„)afc, V = iF54l<'^ + Fs.a,. 

(12) 

With the matrix Riccati equation in Eq. ^ cast in 
terms of Pauli spinors together with coefRcients z^ — 
Z4,Zi: z = 241^^) — iziai, it takes the form 



(13) 



(As an alternative, V and z can also be rendered in terms 
of quaternions (1,— icr^).) and 72 i'^ Eq. (H]) become 



equal and proportional to a unit matrix, (1 



The structure of Eq. admits to the four quantities 
z being real. The effective Hamiltonian in Eq. (jlOp in 
terms of these z becomes 



H^|j'2) ^ h(1'2) _ e,JkZ^F^Jak T F^jZ^Uj ± F^^z.a,. (14) 

We can now construct a five-dimensional unit vector 
rfi out of the four z. 



'2z„ 



(1 + ^^)' 



{I + -4)'' 



H,v=l-A. (15) 



The nonlinear Eq. ([HI), or Eq. in z, becomes of simple, 
linear Bloch-like form, 



IV. THE GENERAL SU(4) HAMILTONIAN 
INVOLVING ALL FIFTEEN OPERATORS 

Instead of the Hamiltonians considered in Section III 
which involve sub-algebras of the full two-qubit system, 
consider an arbitrary 4x4 Hamiltonian with its entire 
complement of 15 operators/matrices. Such a is ob- 
tained by adding to the previous Spin(5) Hamiltonian 

considered above the five additional terms, Fq^u^-P + 
-F640'i^^ + -FejCTa^^CTp^ Correspondingly, Eq. ([TTl) gets an 
additional term zti^gsl'^' in the diagonal H*^^^^) while in 
V, the F^^ are replaced by _F5^ — iFq^. Thus, the full 
SU(4) amounts to a simple modification of the previously 
considered Spin(5) by adding a term proportional to the 
unit operator to the diagonal blocks and making the four 
i^5yj complex, with i^g/^ absorbed as their imaginary parts. 
The Riccati Eq. (fTH]) . now for complex z, becomes 



Zt, = F5^{l-zl)~iFe^{l + zl) + 2F^^z^ 

+ 2(F5^ + iFe^)z„Zf, - 2iF65Z^,/i, 1/ = 1 - 4.(17) 

The two gammas in Eq. ^ are given by 



7i,2 = {1 + zl)l^^^ +i{z*Z4- zlz,)a, 
± -ieijkiztZ* - ZjZ*)ak- 



(18) 



2^ IF^jjj 



/i, 1/ = 1 



(16) 



As in the single spin case, this represents an inverse 
stereographic projection, now from the four-dimensional 
plane z G i?^ to the four-sphere 5"'. It provides a higher- 
dimensional polarization vector for describing such two 
spin problems. With z so described, the two effective 
SU(2) Hamiltonians in Eq. p^ . when solved in turn, 
give the complete solution. In all, such Hamiltonians 
possessing Spin(5) symmetry are, therefore, described by 
the geometrical picture of one S^ and two S^ spheres 
along with two phases. 

su(3) sub-algebra: Four-level systems with only two 
independent energy parameters along the Hamiltonian's 
diagonal and three complex off-diagonal couplings con- 
stitute a su(3) sub-algebra with 8 parameters. A general 
three-level system, embedded into four with the fourth 
level completely uncoupled, constitutes a trivial example 
of such an su(3) sub-algebra but less trivial examples can 
also occur. The z now has two non-zero complex z for 
a total of four parameters. The description of this four- 
dimensional manifold, as well as the remaining SU(2) and 



Their square-root matrices gi^2 can also be evaluated in 
terms of the Pauli matrices and the two SU(2) effective 
Hamiltonians then constructed in explicitly traceless and 
Hermitian form. 

Just as the very structure of Eq. (fT3| suggests that z^ 
and (1 — z^) with suitable normalization define a five- 
dimensional unit vector m in Eq. ([TS]), the occurrence of 



i, (1 — zl), (1 -I- zl) in Eq. p?)) suggests now the intro- 



duction of six quantities according to 



-2z^ 



(1 



with 



me 



*2xl/2 



(19) 



D EE (l + 2|z^|^ + z^zr) 
= -2F65 + iF5^{z;^z^) + Fe^{z; + z,,). (20) 

As with the so(5) case in Section HI, with such a set of 
six complex quantities m, the nonlinear Riccati equation 
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for the four complex in Eq. (|17p becomes a linear 
Bloch-like equation as before, 

■m^^2Ffj_ymy, ^,iy=l-6. (21) 

Once again, the obey a first-order equation with an 
antisymmetric matrix which describes rotations. Since 
the 15 F^i, are real, the real and imaginary parts of 
the six each obey such a rotational transformation. 
These six-dimensional rotations reflect the isomorphism 
between the groups SU(4) and S0(6) (more accurately, 
its covering group Spin(6)) and suggest a mapping be- 
tween their generators (see Appendix B). 

To get a geometrical picture of the manifold to, we 
note first the relations, 



such a manifold by one dimension to [SU(4)/SU(2) x 
SU(2) X U(l)], which is equivalent to the reduction from 
the Stiefel to a Grassmannian manifold G(4, 2, C) ac- 
cording to St(6, 2, R) 9i G(4,2,C)x U(l). Such a Grass- 
mannian manifold, which has eight dimensions, thereby 
describes the z in Eq. Q or its equivalent in Eq. (fT7|) 
or TO^ in Eq. (fT9|) . 

A more accessible geometrical picture is to consider a 
single five-sphere S^ embedded in six-dimensional space 
and two six-dimensional unit vectors from the origin to 
the surface to represent the real and imaginary parts of 
TO^. The two vectors are always taken as orthogonal, so 
that one views such an orthogonally-coupled pair rotat- 
ing within the sphere (28l |. This nine-dimensional object, 
combined with the zero reference of (j) being undefined, 
is our eight-dimensional manifold of interest. 



w'=0, |to^|2 = 2, (22) 

which amount to three constraints. In addition, only the 
derivative, not the value, of (j) is determined in Eq. ((20)) . 
Thereby, the number of independent parameters in 
is eight just as in the complex z^, themselves built from 
z. The description of such an eight-dimensional manifold 
will be taken up in the next sub-section but we note here 
the reduction to the previous so(5) example. This follows 
upon setting Fgs = 0, Fg^ = which makes (p = and 
D = {1 + z1) in Eq. (^0]) . and reduces to^ and tos to 
the values in Eq. (|15p whereas TOg = —i. This, of course, 
makes m a five-dimensional unit vector and its manifold 
the four-sphere S''. The first relation in Eq. (P^ . of the 
vanishing of a square, hints at Grassmannian elements, 
to be discussed further below. 



B. Description in Pliicker coordinates 

An alternative view of these manifolds is pro- 
vided in terms of what are termed Pliicker coordi- 
nates, defined as a set of six complex parameters 
(Pi2, Pi3, Pi4, -P23, -P24, ^34) formed as minors of the 2x4 
sub-matrix of the last two columns of an arbitrary, uni- 
tary SU(4) matrix [13], 



/ Uii U12 Ui3 Ui4 \ 

U21 U22 U23 U24 I 

"31 ""32 U33 U34 I 

\ U41 U42 U43 U44 ) 



A. Nature of the manifold describing {z, m) for 
general SU(4) 

Our construction of the evolution operator for (N ~ 
4,n — 2) in Eq. ((3]) is in terms of the eight-dimensional 
base manifold z and a fiber consisting of two residual 
SU(2) along its diagonal blocks and a U(l) phase between 
them: SU(4) ^ [SU(4)/SU(2) x SU(2) x U(l)] x [SU(2) 
X SU(2) x U(l)]. To describe the former base manifold, 
consider first [SU(4)/SU(2) x SU(2)], which is a nine- 
dimensional manifold. It can also be described in terms 
of spin-groups as Spin(6)/Spin(4). The six complex 
in Eq. (fT8|) with the three constraints in Eq. ([22]) consti- 
tute such a manifold called a Stiefel manifold St(6, 2, R) 
= this name being given to manifolds consisting of n 
orthogonal vectors from an iV-dimensional space [131 ■ 
Geometrically, the second relation in Eq. ([22]) states that 
the real and imaginary parts of to are six-dimensional 
unit vectors while the first relation expresses their mutual 
orthogonality. Therefore, one can view the manifold as 
a five-sphere S^ with another four-sphere S^ attached at 
each point on it. The absolute value of the phase param- 
eter (j) in Eq. (|19p and Eq. (|20p being undefined, reduces 



They obey the relations 



P12P34 - P13P24 + P14P23 = 0, ^ \P^J P = 1. (24) 
They are combinations of the to^ according to 



/ P12 






/ 


ITOg - 


- m5 


\ 


Pis 








imi - 


- TO2 




Pl4 




1 




—im^ 


+ m4 




P23 




" 2 




—im^ 


— 7714 




P24 








— ITOl 


+ TO2 




\P34 






I 




- 


) 



(25) 



The linear equations for to^ in Eq. (|20p translate into 
a similar linear equation 



ZP = HP, P = (Pi2,-Pl3,Pl4,P23,P24,P34), (26) 

with 
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H 



/ ^^11, 22 


H41 


H31 


—H42 


H32 









^11,33 




—H12 





^H^2 


Hl3 


— H43 


^^11,44 





H12 


H42 




— H24 


— H21 





H22.33 


H34 











H21 




H22M 


—H41 




\ 


— H23 


H24 


— iJl3 


— Hi4 


H33AA 


) 



p — 



(27) 

where we have adopted the notation for the diagonal en- 
tries: Hii^jj = Hii + Hjj. 

Actually, the above equations for P can be arrived at 
directly from the evolution equation iU = HU because 
the elements of P are quadratic in the elements of U in 
Eq. (El: P, 



z can be defined in terms of the two minors on the right 
in Eq 



and iiif^ = Hkjuf \ Also, 



"13 Ml4 
U23 U24 



/ 



M33 M34 
M43 M44 



(28) 



the matrix in the denominator assumed to be non- 
singular. Writing U in Eq. ((23)) in the form in Eq. ([3]), 
the first factor Ui involving z is a map of the Grass- 
mannian manifold G(4, 2, C) onto C^, and provides 
a partial coordinization of that manifold. Elements of 
G(4, 2, C) are two-dimensional complex hyperplancs 
spanned by vectors U3 = (wi3, U23, M33, "43)"^ and U4 = 
(mi4, U24, W34, "44)^. The Pliicker coordinates provide a 
unique identification of such planes. They are an analog 
of the coordinization of the n-dimensional sphere S" by 
an [n + l)-dimcnsional unit vector rfi as in Section III. 

The matrix Hp in Eq. (|77|) being Hermitian, P^P 
= constant = 1. This can be verified by the relation 
between P's and m's in Eq. (P5|) which involves a uni- 
tary matrix so that P^P — ^m^m, and combining with 
Eq. ([221) • Further, a symplectic structure can be intro- 
duced. Defining a 6 x 6 matrix = with non-zero 
entries of 1 only along the anti-diagonal, the first rela- 
tion in Eq. ([24]) can be rendered as P-^JIP = 0, and the 
matrix Hp, a generator of the symplectic group Sp(6, 
C), 



HpO + rtii^p = Tr(Hp)rj = 0. 



(29) 



Any two vectors P;, evolving according to Eq. ([26|l . satisfy 
Pf (i)nP2(i) ^constant. If PiJ7P2 = 0, then the two 
hyperplanes defined by P; intersect, and if |Pir2P2| = 1, 
they do not. 

Geometrically, the set P^^P = 1 is a sphere S^^, the alge- 
braic relation PJIP = determining a 9-dimensional sub- 
manifold, an intersection between S^^ and the affine variety 
of roots of the polynomial equation PfJP = 0. This mani- 
fold may be denoted Multiplication by a phase acts as 
a transformation group on this manifold, that is, if P e 5R^, 
then Ve^'t' e Therefore, G(4, 2, C) is a quotient 

space K^/U(l) and has eight dimensions. The connection 



to SU(4) is, as noted before, 3?^ ?^SU(4)/(SU(2) x SU(2)) 
= Spin(6)/ Spin(4). The stability sub-group of a vector 
P e is SU(2) X SU(2) while the stability sub-group of 
3?7U(1) isSU(2) X SU(2) X U(l). Since Spin(6)/Spin(5) 
= and Spin(5)/Spin(4) = S^, we can identify the fibra- 
tion of with S^x S^. 



V. SUMMARY 

We have presented a complete analysis of the evolution 
operator for SU(Af), setting up its construction in a hier- 
archical way in terms of those for smaller SU(A^ — n) and 
SU(n), with n < N and arbitrary. The evolution operator 
is written as a product of two N x N matrices, the second 
of which is block diagonal in {N — n) x (N — n) and n x n 
of the smaller groups. The first factor is obtained through 
a z, which is an {N — n) x n complex matrix obeying a ma- 
trix Riccati equation. Its solutions determine both the first 
factor as well as the Hermitian matrices for the subsequent 
N — n and n evolution problems. 

This general constructive method is applied especially to 
a four-level system with special emphasis on two qubits. 
The general symmetry is of SU(4), a 15-parameter group. 
Our procedure expresses the evolution operator as a prod- 
uct of two 4x4 matrices, the second of which is block 
diagonal, each block an SU(2) problem. The z is also a 
2x2 matrix with complex entries in general and obeys a 
matrix Riccati equation. Alternatively, we transform z into 
a six-dimensional complex vector to, whose real and imag- 
inary parts both separately undergo linear, six-dimensional 
rotational transformations. This is exactly analogous to the 
linear Bloch equation for real three-dimensional rotations of 
a vector to represent the evolution operator for a single spin 
in a magnetic field. 

Just as a Bloch sphere describes the three-dimensional 
vector TO for a single spin (and, together with a phase, the 
complete SU(2)), we also present the geometrical manifold 
describing z or its equivalent six-dimensional complex vector 
TO. Together with two residual SU(2) problems and a phase, 
this provides a complete description of the quantum evolu- 
tion operator for SU(4). For certain sub-algebras of SU(4), 
the manifold is an analogous higher-dimensional sphere; a 
four-sphere, for example, for an so(5) sub-algebra. For the 
most general SU(4), we have an eight-dimensional Grass- 
mannian manifold. We provide a picture of it as two five- 
spheres with an orthogonality and phase constraint. These 
geometrical objects may serve for all possible four-level and 
two qubit systems the useful purpose that the Bloch sphere 
has for two-level and single qubit problems in physics. 



APPENDIX A: EXTENSION TO NON-UNITARY 
EVOLUTION FOR A NON-HERMITIAN 
HAMILTONIAN 

The iterative method of Section II for the evolution op- 
erator in Eq. ([2]) through writing it as in Eq. ([3]) applies 
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also when H in Eq. ([l} is not Hermitian and, therefore, the 
evolution not unitary. However, z and w in Eq. Q are no 
longer simply related as in Eq. (|4]) but obey independent 
equations, the former still in Riccati form but the latter 
given in terms of z. Thus, instead of Eq. ([1]), consider 



H(^-")(t) Y{t) 
Yt(t) H(«)(<) 



(30) 



where we have again indicated by tildes non-Hermiticity, 
and V and Y are not equal but independent. 

Writing UW(t) again as in Eq. ([3]), Eq. ([6]) now becomes 

ii = (H(^-i^z-zH(") -zYtz + V, 
iW^ = wt(zYLH(^-i))+ (h(")- Ytz)wt+Yt.(31) 

The residual problems of [N — n) and n dimension then 
become 



which is explicitly anti-symmetric. We thus have H = 
2F^^Liyf^. Analogous to the familiar triplet of angular mo- 
mentum generators, six-dimensional generators of S0(6) are 
given by L^^ = —ilf_iu, where the I are 15 real antisymmet- 
ric 6 x 6 matrices with only two non-zero entries, -f 1 in the 
{liiy) and —1 in the {lyfi) position: 



(34) 



Their commutators close: 



■U(N-n) 



H(^-i) - zYt , 
Ot H(") + Ytz ' 

(32) 



APPENDIX B: DESCRIPTION OF EVOLUTION 
AS SIX-DIMENSIONAL ROTATIONS 



For a single spin or qubit, the rewriting of the quantum 
evolution operator, which is complex, as rotational trans- 
formations of a real, unit vector in three dimensions given 
by the Bloch equation, rests on the isomorphism of the 
group SU(2) to S0(3) (or its double covering Spin(3)). A 
similar isomorphism between the groups SU(4) and S0(6) 
(or its extension Spin(6)) underlies the construction in Sec- 
tions III and IV of the complex evolution operator for two 
qubits in terms of rotations of a vector in six dimensions. 
Both groups are described by 15 real parameters through 
an antisymmetric F^^, /i, = 1, 2, ... 6. In Sections III and 
IV, explicit expressions are given for the Hamiltonian with 
each of these parameters multiplying one of the 15 complex 
generators of SU(4) in a standard representation of Pauli 
matrices, (t'^) (XiXf^',!^^) (g)!?'^)^ ^(i) 0^(2) ^^1 alternative 

rendering in terms of the 15 generators of S0(6) is useful 
and recorded here. 

The Hamiltonian in Section IV, apart from a factor of i, 
can be cast in terms of a matrix array 

ai'^ai'^ ai'^ai'^ -^^^ ai'^ 
-ai^Vi^) -a^^a^^ -ai^Vi^) 

(33) 



so that form an so(6) algebra. 

The array in Eq. (|33p is also a convenient display of the 
generators of the various sub-groups of S0(6) of lower- 
dimensional rotations. Either upper left or lower right cor- 
ner 2x2 blocks describe the S0(2) generator of one of the 
qubits. Adding a third row and column gives the full triplet 
of S0(3) generators. To this can be added a next row and 
column of three non-zero entries to give the six generators 
of S0(4). For this purpose, any of the three remaining 
row/column can be employed, each giving an S0(4), the 
three added entries transforming as a vector under S0(3). 
This continues. Adding another row and column's four new 
entries, which transform as a vector under S0(4) (further 
subdividing into three components that transform as a vec- 
tor and one as a scalar under the previous S0(3)), gives the 
ten S0(5) generators. The final sixth row/column adds five 
entries, an S0(5) vector, to give the full 15 generators of 
S0(6). This hierarchical nesting of SO sub-groups, together 
with the corresponding Clifford structure with Pauli matri- 
ces in Eq. p3p . accounts for the richness of the structures 
in the isomorphic groups SU(4) and S0(6), one we have ex- 
ploited in Sections III and IV. Note that the linear Bloch-like 
: equation for to in Eq. (|2ip for a general SU(4) Hamiltonian 
reduces to the same antisymmetric form for its sub-groups 
such as in Eq. psp . all the way down to the standard Bloch 
equation for a single qubit, whose S0(3) antisymmetric 
is usually written as a vector product with a magnetic field. 



9 



[*] Email: arau@phys.lsu.edu 

[1] W. P. Schleich, Quantum Optics in Phase Space (Wiley- 
VCH, Berlin, 2001). 

[2] A. Abragam, The Principles of Nuclear Magnetism 
(Clarendon, Oxford, 1978); R. R. Ernst, G. Bodenhausen, 
and A. Wokaun, Principles of Nuclear Magnetic Reso- 
nance in One and Two Dimensions{C\arendon , Oxford, 
1987). 

[3] H. F. Jones, Groups, Representations and Physics (lOP 
Publishing, Bristol, 1990); J. J. Sakurai, Modern Quantum 
Mechanics (Addison- Wesley, Reading, MA, 1994), Chapter 
3; U. Fano and A. R. P. Rau, Symmetries m Quantum 
Physics (Academic Press, San Diego, 1996). 

[4] K. Blum, Density matrix theory and applications (Plenum, 
New York, 1996); J. J. Sakurai, Modem Quantum Mechan- 
ics (Addison-Wesley, Reading, MA, 1994), Sec. 3.4. 

[5] Kurt Gottfried and Tung-Mow Yan, Quantum Mechan- 
ics: Fundamentals (Springer, New York, 2004), second ed., 
p. 447-8. 

[6] M. Nielsen and I. Chuang, Quantum Computation and 
Quantum Information (Cambridge Univ. Press, New York, 
2000). 

[7] E. Gerjuoy, Phys. Rev. A 67, 052308 (2003); Z. Sun, X. 

Wang, and C. P. Sun, ibid, 75, 062312 (2007); F. Lastra, 

G. Romero, C. E. Lopez, M. F. Santos, and J. C. Retamal, 

ibid, 75, 062324 (2007). 
[8] C. Bemporad, G. Gratta, and P. Vogel, Rev. Mod. Phys. 

74, 297 (2002); M. C. Gonzalez-Garcia and Y. Nir, ibid75, 

345 (2003). 

[9] P. Guttinger, Z. Phys. 73, 169 (1931); E. Majorana, 
Nuovo Cimento 9, 43 (1932); J. Schwinger, Phys. Rev. 51, 
648 (1937); J. J. Sakurai, Modem Quantum Mechanics 
(Addison-Wesley, Reading, MA, 1994), Sec. 3.8. 

[10] I. M. Gelfand, M. M. Kapranov, and A. V. Zelevinsky, Dis- 
criminants, Resultants and Multidimensional Determi- 
nants (Birkhauser, Boston, 1994), chapter 3; Peter Levay, 
J. Phys. A: Math. Gen. 38, 9075 (2005); Phys. Rev. A 71, 
012334 (2006); Phys. Rev. D 74, 024030 (2006). 

[11] A. R. P. Rau, Phys. Rev. Lett. 81, 4785 (1998). 

[12] D. B. Uskov and A. R. P. Rau, arXiv:quant-ph/0511192 and 
Phys. Rev. A 74, 030304 (R) (2006); see also S. Chaturvedi, 
E. Ercolesi, G. Marmo, G. Morandi, N. Mukunda, and R. 
Simon, Pramana (2007). 



[13] A. R. P. Rau and D. B. Uskov, |arXiv:quant-ph/0512135| and 

Phys. Scr. 74, C31 (2006). 
[14] T. W. Reid, Riccati Differential Equations, Mathematics 

in Science and Engineering, Vol. 86 (Academic, New York, 

1972). 

[15] S. Lang, Fundamentals of Differential Geometry, Gradu- 
ate Texts in Mathematics, Vol. 191 (Springer, New York, 
1999). 

[16] See, for instance, F. E. Close, An Introduction to Quarks 

and Partons (Academic, London, 1979). 
[17] F. J. Ohkawa, J. Phys. Soc. Japan 52, 4785 (1998). 
[18] A. R. P. Rau, Phys. Rev. A 61, 032301 (2000). 
[19] F. J. M. van de Ven and C. W. Hilbers, J. Magn. Reson. 

54, 512 (1983). 

[20] J. Novotny, G. Alber, and I. Jex, Phys. Rev. A 73, 062311 
(2006). 

[21] A. R. P. Rau, G. Selvaraj, and D. Uskov, Phys. Rev. A 71, 
062316 (2005). 

[22] T. Yamamoto, Yu. A. Pashkin, 0. Astafiev, Y. Nakamura, 
and J. S. Tsai, Nature 425, 941 (2003). 

[23] F. Legare, Phys. Rev. A 68, 063403 (2003); V. Delgado 
and J. M. Gomez Llorente, ibid 68, 022503 (2003). 

[24] The SO(n) orthogonal group of rotations is familiar and 
suffices for most of our purposes, although technically it 
is the Spin(n) groups which are the double covering groups 
that occur as sub-groups of SU(n) with the Clifford structure 
of the Pauli matrices. Even more accurately, in the sense of 
the Cartan decomposition, it is SU(2) ® SU{2)/Z2, where 
Z2 is a group of ±1, which is isomorphic to S0(4) whereas 
SU(2) ® SU(2) is isomorphic to Spin(4). 

[25] R. Cabrera, C. Rangan, and W. E. Baylis, Phys. Rev. A 76, 
033401 (2007). 

[26] Sai Vinjanampathy and A. R. P. Rau, unpublished (2007). 

[27] Ingemar Bengtsson and Karol Zyczkowski, Geometry of 
Quantum States (Cambridge Univ. Press, 2006). 

[28] A lower dimensional analog is instructive. Two real quanti- 
ties, {xi,X2) with xf-'rxi ~ 1 define a circle or one-sphere 
S^. Two complex numbers (2:1,22), with z — x-\-iy, obey- 
ing the constraints in Eq. (|22p. \zi\'^ = 1, zf = 0, also 
define a one-dimensional manifold of two real, orthogonal, 
unit vectors [x,y). The in Eq. (|19p and Eq. I|22p are 
the six-dimensional extension of this. 



